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This  paper  investigates  the  acoustic  uncertainty  due  to  hydrographic  data  error  and  in 
turn  to  determine  the  necessity  of  a  near  real  time  ocean  analysis  capability  such  as  the 
Naval  Oceanographic  Office’s  (NAVOCEANO)  Modular  Ocean  Data  Assimilation 
System  (MODAS)  model  in  shallow  water  (such  as  the  Yellow  Sea)  mine  hunting 
applications  using  the  Navy’s  Comprehensive  Acoustic  Simulation  System/Gaussian 
Ray  Bundle  (CASS/GRAB)  model.  To  simulate  hydrographic  data  uncertainty, 
Gausian-type  errors  (produced  using  the  random  number  generator  in  MATLAB)  with 
zero  mean  and  three  standard  deviations  (1  m/s,  5  m/s,  and  10  m/s)  are  added  to  the 
sound  profile.  It  is  found  that  the  acoustic  uncertainty  depends  on  the  location  of  the 
error  and  sound  sources.  It  is  more  sensitive  to  errors  in  the  isothermal  structure  in  the 
winter  than  in  the  layered  structure  in  the  summer. 

1  Introduction 

The  major  threats  in  the  littoral  are  diesel  submarines  and  sea  mines.  The 
combination  of  improvements  in  noise  reducing  technology  and  the  development  of 
Air  Independent  Propulsion  (AIP)  technology  have  made  diesel  submarines  very 
difficult  to  detect  in  both  the  littoral  and  blue  waters.  After  a  weapon  platform  has 
detected  its  targets,  the  sensors  on  torpedoes  designed  for  blue  water  operations  are 
not  designed  to  acquire  a  target  in  a  reverberation-crippling  environment.  Recently, 
the  U.S.  Navy  has  focused  much  of  its  research  and  development  efforts  in 
designing  high  frequency  sensors  and  corresponding  acoustic  models  to  overcome 
the  threat  in  the  littoral.  The  Comprehensive  Acoustic  Simulation  System  (CASS) 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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using  the  Gaussian  Ray  Bundle  (GRAB)  model  is  a  valuable  tool  for  the  AN/SQQ-32 
mine  hunting  detection  and  classification  sonar.  The  performance  of  this  model,  as  in  all 
models,  is  determined  by  the  accuracy  of  its  inputs  such  as  sea  surface  conditions, 
bathymetry,  bottom  type,  and  sound  speed  profiles.  Here,  the  effect  of  sound  speed  errors 
(i.e.,  hydrographic  errors)  on  the  acoustic  uncertainty  in  the  Yellow  Sea  is  investigated  using 
CASS/GRAB. 

2  Environment  of  the  Yellow  Sea 


2. 1  Geology  and  Structure 

The  Yellow  Sea  is  a  semi-enclosed  basin  situated  between  China  and  the  Korean 
peninsula  with  the  Bohai  Sea  to  the  northwest  and  the  East  China  Sea  to  the  south.  The 
Yellow  Sea  is  a  large  shallow  water  basin  covering  an  area  of  approximately  295,000 
km2.  The  water  depth  over  most  of  the  area  is  less  than  50  m  (Fig.  1).  The  bottom 
sediment  of  the  central  and  western  regions  of  the  Yellow  Sea  consists  primarily  of  mud 
and  the  eastern  region  is  primarily  sand.  The  mud  sedimentation  in  the  central  and 
northwestern  regions  of  the  Yellow  Sea  is  due  to  the  runoff  from  the  great  rivers  of 
China.  Four  regions  with  different  bottom  types  were  selected  for  the  acoustic  model 
runs  in  this  study  (Fig.  2):  (a)  rock  bottom  type  which  is  located  in  the  north-central 
Yellow  Sea  at  37°-37.5°N,  123°-123.8°  E,  (b)  gravel  bottom  type  which  is  located  in  the 
northern  Yellow  Sea  at  38.4°-39°N,  122°-123°E,  (c)  sand  bottom  type  which  is  located 
in  the  southeastern  Yellow  Sea  at  35.5°-36.5°  N,  124.5°-126.2°  E,  and  (d)  mud  bottom 
type  which  is  located  in  the  south-central  Yellow  Sea  at  35°-36.5°N,  123°-  124.5°  E. 


Fig.  1.  Bottom  Topography  of  the  Yellow  Sea 

and  the  surrounding  regions  Fig.  2.  Yellow  Sea  Bottom  sediment 

chart 


2.2.  Oceanography 

The  four  seasons  in  the  Yellow  Sea  are  defined  as  follows:  the  winter  months  run  from 
January  through  March;  the  spring  months  run  from  April  through  June;  the  summer 
months  run  from  July  through  September;  and  the  fall  months  run  from  October  through 
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December.  The  two  main  characteristic  temperature  profiles  of  the  Yellow  Sea  are 
during  the  winter  and  the  summer  months.  In  the  winter  months,  the  temperature 
profiles  throughout  the  region  are  characterized  as  isothermal  (Fig.  3a).  In  the  summer 
months,  the  temperature  profiles  throughout  the  region  are  characterized  by  a  multi-layer 
profile  consisting  of  a  mixed  layer,  a  thermocline,  and  a  deep  layer  (Fig.  3b). 


Fig.  3.  Eastern  Yellow  Sea  (around  36  N)  temperature  profiles  during  1950-1988;  (a)  January  and  (b) 
June.  Solid  dots  show  the  location  of  the  observation  stations  (From  Chu  et  al.  1997b). 


3  CASS/GRAB  Model 

CASS/GRAB  is  an  active  and  passive  range  dependent  propagation,  reverberation,  and 
signal  excess  acoustic  model  and  is  accepted  as  the  Navy’s  standard  model.  The  GRAB 
model’s  main  function  is  to  calculate  eigenrays  in  range-dependent  environments  in  the 
frequency  band  600  Hz  to  100  kHz  and  to  use  the  eigenrays  to  calculate  propagation 
loss.  The  CASS  model  is  the  range  dependent  improvement  of  the  Generic  Sonar  model 
(GSM).  CASS  performs  range  independent  monostatic  and  bistatic  active  signal  excess 
calculations.  The  major  difference  between  the  GRAB  model  and  a  classic  ray  path  is 
that  the  amplitude  of  the  Gaussian  ray  bundles  is  global,  affecting  all  depths  to  some 
degree,  whereas  classic  ray  path  amplitudes  are  local.  GRAB  calculates  amplitude 
globally  by  distributing  the  amplitudes  according  to  the  Gaussian  equation 


Gv=  (0.5)(max(Az,47T?L)), 

where  the  Ty  represents  losses  due  to  volume  attenuation  and  boundary  interaction;  gv 

is  the  effective  standard  deviation  of  the  Gaussian  width;  and  jBv <0  is  a  factor  that 
depends  only  on  the  source  and  is  chosen  so  that  the  energy  within  a  geometric-acoustic 
ray  tube  equals  the  energy  within  a  Gaussian  ray  bundle.  The  variable  zv  is  the  depth 
along  the  vth  test  ray  at  range  r,  z  is  the  target  depth,  pr  is  the  horizontal  slowness,  Az  is 
the  change  in  ray  depth  at  constant  range  due  to  a  change  in  source  angle,  and  X  is  the 
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wavelength.  The  selection  of  the  effective  standard  deviation  gv  is  the  weakest 
component  in  providing  a  firm  theoretical  basis  for  the  GRAB  model.  GRAB  computes 
the  random  or  coherent  propagation  loss  from  the  eigenrays  stored  in  the  eigenray  file 
and  stores  in  them  in  separate  pressure  files  (Aidala  et  al.  1998). 


4  Modular  Ocean  Data  Assimilation  System 

Modular  Ocean  Data  Assimilation  System  (MODAS),  recently  developed  at  the  Naval 
Research  Laboratory  (NRL),  uses  a  modular  approach  to  generate  three-dimensional 
gridded  fields  of  temperature  and  salinity.  Its  data  assimilation  capabilities  may  be 
applied  to  a  wide  range  of  input  data,  including  randomly  located  in-situ,  satellite,  and 
climatological  data.  Available  measurements  from  any  or  all  of  these  sources  are 
incorporated  into  a  three-dimensional,  smoothly  gridded  output  field  of  temperature  and 
salinity.  MODAS’  primary  outputs  are  temperature  and  salinity  fields  that  may  be  used 
to  calculate  three-dimensional  sound  speed  fields.  The  sound  speed  field,  in  turn,  may 
be  used  to  drive  acoustic  performance  prediction  scenarios,  including  simulations, 
tactical  decision  aids,  and  other  capabilities.  These  are  employed  in  a  wide  variety  of 
naval  applications  and  tactical  decision  aids  (TDAs)  (Fox  et  al.  2002). 


5  Acoustic  Characteristics 

On  February  15,  2000,  the  sound  speed  profile  at  36.4°N,  124.4°E  (mud  bottom) 
from  MODAS  is  quite  uniform.  The  sound  speed  decreases  0.2  m/s  from  the  surface  to 
8.2  ft  depth,  increases  0.1  m/s  from  8.2  ft  to  41.0  ft  depth,  and  0.2  m/s  from  41.0  ft  to 
57.4  ft  depth.  The  weak  sound  speed  minimum  at  near  surface  (8.2  ft  depth)  generates  a 
weak  sound  channel  with  sound  source  at  25  ft  (Figure  4)  and  125  ft  (Figure  5). 


Table  1.  Sound  speed  from  MODAS  on  February  15,  2000,  at  36.4°  N,  124.4°  E, 

mud  bottom 

Depth  (Feet) 

M/S 

0.00 

1479.90 

8.20 

1479.70 

24.60 

1479.80 

41.00 

1479.80 

57.40 

1480.00 

82.00 

1480.00 

106.60 

1480.10 

131.20 

1480.30 

164.00 

1480.50 

205.00 

1480.40 

246.00 

1480.40 

6  Acoustic  Uncertainties 

To  simulate  hydrographic  data  uncertainty,  a  Gausian-type  error  (produced  using  the 
random  number  generator  in  MATLAB)  with  zero  mean  and  three  standard  deviations 
(1  m/s,  5  m/s,  and  10  m/s)  is  added  to  the  MODAS  sound  profiles.  All  the  sound  speed 
profiles  with  the  mud  bottom  are  selected.  Two  sets  of  hydrographic  data  (MODAS  and 
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MOD  AS  with  error)  are  inputted  into  the  CASS/GRAB  model.  The  model  was 
integrated  with  two  sound  source  depths  (25  and  125  ft)  and  two  seasons  (February  15, 
2000  and  August  15,  2000)  to  capture  the  effect  of  the  hydrographic  error  on  the 
acoustic  uncertainty  in  the  Yellow  Sea. 

Histograms  of  the  detection  range  difference  (acoustic  uncertainty)  between  control 
(no  error)  and  sensitivity  (with  error)  runs  demonstrate  that  the  acoustic  uncertainty  has 
non-Gaussian-type  distribution  in  winter  (Figure  6)  and  Gaussian-type  distribution  in 
summer  (Figure  7)  and  is  much  larger  in  winter  than  in  summer.  This  indicates  that  the 
isothermal  structure  of  the  winter  profiles  is  much  more  susceptible  to  errors  in  sound 
speed.  The  acoustic  uncertainty  depends  on  where  the  random  error  is  situated  in  the 
water  column  in  relation  to  the  position  of  the  source.  For  a  specific  profile,  if  an  error 
of  1  m/s  is  positioned  within  approximately  5  feet  of  the  source  depth  and  an  error  of  10 
m/s  is  positioned  greater  than  the  5  feet  of  the  source  depth,  the  1  m/s  error  will  have  a 
much  greater  effect  on  the  acoustic  transmission.  If  the  error  near  the  source  is  positive, 
the  gradient  that  is  formed  in  the  sound  speed  profile  will  decrease  detection  ranges. 
If  the  error  is  negative,  the  gradient  that  is  formed  in  the  sound  speed  profile  will 
increase  detection  ranges. 


Figure  4.  Acoustic  transmission  with  source  depth  of  25  ft  on  February  15,  2000  at  36.4  N  124.4  E  (mud 
bottom)  using  the  MOD  AS  data:  (a)  ray  trace,  and  (b)  signal  excess.  Notice  that  a  maximum  detection 
range  is  260  yd  near  the  sound  source. 


Figure  5.  Acoustic  transmission  with  source  depth  of  125  ft  on  February  15,  2000  at  36.4  N  124.4  E 
(mud  bottom)  using  the  MOD  AS  data:  (a)  ray  trace,  and  (b)  signal  excess.  Notice  that  a  maximum 
detection  range  is  145  yd  near  the  sound  source. 
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Figure  6.  Histograms  of  the  detection  range  difference  caused  by  the  Gaussion-type  errors  in  the  sound 
speed  profiles  on  February  15,  2000  with  the  mud  bottom  and  125  ft  source  depth:  (a)  MODAS  minus 
MOD  AS  with  1  m/s  error,  (b)  MODAS  minus  MODAS  with  5  m/s  error,  and  (c)  MODAS  minus 
MODAS  with  10  m/s  error. 


Figure  7.  Histograms  of  the  detection  range  difference  caused  by  the  Gaussion-type  errors  in  the  sound 
speed  profiles  on  August  15,  2000  with  the  mud  bottom  and  125  ft  source  depth:  (a)  MODAS  minus 
MODAS  with  1  m/s  error,  (b)  MODAS  minus  MODAS  with  5  m/s  error,  and  (c)  MODAS  minus 
MODAS  with  10  m/s  error. 

7  Hydrographic  Errors  at  the  Source  Depth 


When  an  error  (+1  m/s)  was  added  into  the  MODAS  sound  speed  profile  at  both  source 
depths  (25  and  125  ft),  a  shadow  zone  was  formed  in  front  of  the  source  that 
significantly  decreased  the  detection  ranges  at  that  depth  (Figures  8  and  9).  When  an 
error  (-lm/s)  was  added  into  the  MODAS  sound  speed  profile  at  both  source  depths,  a 
strong  sound  channel  formed  that  dramatically  increased  detection  ranges  at  that  depth 
(Figures  10-11). 


Figure  8.  MODAS  with  +1  m/s  sound  speed  error  at  the  source  depth  (25  ft)  on  February  15,  2000,  36.4 
N  124.4  E,  mud  bottom:  (a)  ray  trace,  (b)  signal  excess  (maximum  detection  range  at  source  depth  =  175 
yd,  =max  detection  range  at  the  source  depth  =  -  85  yd). 
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Figure  9.  MOD  AS  with  +1  m/s  sound  speed  error  at  source  depth  (125  ft)  on  February  15,  2000,  36.4  N 
124.4  E,  mud  bottom:  (a)  ray  trace,  and  (b)  signal  excess  (maximum  detection  range  near  the  source 
depth  =150  yd,  =miax  detection  range  near  the  source  depth  =  -5  yd). 


Figure  10.  MODAS  with  -1  m/s  sound  speed  error  at  the  source  depth  (25  ft)  on  February  15,  2000,  36.4 
N  124.4  E,  mud  bottom:  (a)  ray  trace,  and  (b)  signal  excess  (maximum  detection  range  near  the  source 
depth  >1000  yd,  =miax  detection  range  near  the  source  depth  >740  yd). 


Figure  11.  MODAS  with  -1  m/s  sound  speed  error  at  the  source  depth  (125  ft)  on  February  15,  2000, 

36.4  N  124.4  E,  mud  bottom:  (a)  ray  trace,  and  (b)  signal  excess  (maximum  detection  range  near  the 
source  depth  >1000  yd,  =miax  detection  range  near  the  source  depth  >855  yd). 

8  Conclusions 

(1)  The  seasonal  variation  in  acoustic  transmission  in  the  Yellow  Sea  for  all  regions  was 
mainly  due  to  the  isothermal  structure  in  the  winter  and  a  multi-layer  thermal  structure 
in  the  summer.  The  acoustic  transmission  in  the  winter  is  shorter  due  to  the  effect  of  the 
isothermal  structure  of  the  sound  speed  profile.  The  acoustic  transmission  in  the 
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summer  is  significantly  longer  due  to  the  down  bending  effects  of  the  multi-layer 
structure  of  the  sound  speed  profiles,  which  produce  convergence  zone  and  caustics. 

(2)  The  acoustic  uncertainty  due  to  error  in  sound  speed  profiles  depends  on  location  of 
that  error  relative  to  the  sound  source.  It  has  non-Gaussian-type  distribution  in  winter 
and  Gaussian-type  distribution  in  summer  and  is  much  larger  in  winter  than  in  summer. 

(3)  In  winter,  when  an  error  (+1  m/s)  was  added  into  the  sound  speed  profile  at  the 
source  depth,  a  shadow  zone  was  formed  in  front  of  the  source  that  significantly 
decreased  the  detection  ranges  at  that  depth.  When  an  error  (-lm/s)  was  added  into  the 
sound  speed  profile  at  the  source  depth,  a  strong  sound  channel  formed  that  dramatically 
increased  detection  ranges  at  that  depth. 
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